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Abstract. We study the Fourier transform of polynomials in an orthogonal family, taken with 
respect to the orthogonality measure. Mastering the asymptotic properties of these transforms, that 
we call Fourier-Bessel functions, in the argument, the order, and in certain combinations of the 
two is required to solve a number of problems arising in quantum mechanics. We present known 
results, new approaches and open conjectures, hoping to justify our belief that the importance of 
these investigations extends beyond the application just mentioned, and may involve interesting 
discoveries. 
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1. Introduction and examples. Let /x be a positive measure, for which the 
moment problem is determined, and let {p„(/i; s)}„gN be its orthogonal polynomials. 
The Fourier-Bessel functions (F-B. for short) J7n(/^;i) are the Fourier transforms of 
p„(/i; s) with respect to /i: 

(1.1) Jn{ti;t) := j dfi{s) p„{fi;s) e~''\ 

This nomenclature follows — for lack of better candidates — from the simple observation 
that when /z is the continuous measure with density dfi{s) = ^^^-^==5, and therefore 
PnifJ-] s) are the (properly normalized) Chebyshev polynomials, the F-B. functions are 
the usual integer order Bessel functions: J^nifJ-', t) = {—i)"Jn{t). When the measure is 
symmetrical with respect to the origin, as in this case, the F-B. functions are either 
real, or purely imaginary. A graph of the first few F-B. functions, multiplied by 
z", is displayed in Fig. 11.11 for a singular continuous measure supported on a real 
Julia set (to be introduced in the following). Notice the joyful oscillations that these 
F-B. functions feature, as opposed to the more disciplined, and in the end boring 
attitude of the Jn's. This paper wants to be an ode to the fascinating properties of 
F-B. functions of singular continuous measures, that in my opinion are still largely 
unexplored: I shall present a few results, but mostly open problems. The style of 
this paper will be suggestive of possible developments, rather than assertive of formal 
results, and at times I shall gladly renounce to rigor in favor of intuition, hoping with 
confidence that others will take up where I have left, and complete the picture. In 
this way, I believe to be correctly interpreting Ed's attitude towards mathematics as 
a communal endeavor, and it is not only a pleasure for me, but an honor, to dedicate 
to him these notes. 

The asymptotic of F-B. functions for large values of the argument, t, is a classical 
theme of investigation, especially when n — 0, since ^o{ii;t) is the Fourier transform 
of the measure ^ 012 El S El • In this study, the nature of the orthogonality measure 
H plays a major role. In fact, it is in the realm of singular, multi-fractal measures 

* Dedicated to Ed Saff for his Sixtieth. 

^Center for Non-linear and Complex Systems, Universita dell'Insubria, Via Valleggio 11, 22100 
Como, Italy (giorgioauninsubria.it). 



2 



G. MANTICA 




that the most interesting phenomena appear. First of aU, at difference with the usual 
Bessel case, convergence of J7n(/i;t) to zero is not to be expected, and indeed in 
Figure [L^ that depicts a much larger argument range than Fig. If .11 this time for 
a measure associated with a linear Iterated Function System, bursts of "activity" of 
j7o(/i;t) are observed, amidst zones of quiescence. Because of similarities with the 
theory of turbulence, I have termed this phenomenon and its consequences quantum 
intermittency [HIEIIH]- 

A common technique to cope with these bursts is to take suitable time averages, 
like Cesaro's. After averaging, decay of JnilJ-', t) to zero actually takes place, according 
to an algebraic law. Now, two main problems can be investigated: the decay of the 
averaged F-B. functions themselves, and that of their (averaged) square moduli, this 
second problem having received larger attention than the first. In two recent papers 
pill(J| we have collected known and new results on these questions, under the unifying 
theme of Mellin transforms. The following scheme is encountered in these theorems, 
under very broad hypotheses (typically, the existence of orthogonal polynomials): for 
any x less than the divergence abscissa of a potential theoretic function, the Cesaro 
average of F-B. functions (or of their square moduli) decays faster than t^^. The 
divergence abscissas entering these theorems are identified as the local dimension of 
the measure at zero in the first case, and as the correlation dimension of the measure 
in the second. The appearance of dimensional quantities of the orthogonality measure 
is not accidental: indeed, they play a major role in the asymptotics of F-B. functions, 
as it will become apparent in the following. 

Quite different is the asymptotic behavior for large values of the order, and fixed 
argument. A general result can be obtained on the basis of a Chebyshev expansion of 
the matrix exponential 11 : this theorem states that under the sole hypothesis that 
the support of fi is bounded, at fixed time t, for any a > 0, there exist a constant Ca 
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Fig. 1.2. F-B. function JoCa*;*) for an I.F.S. measure, over a larger t-scale than in Fig. 11. Ip . 

SO that the F-B. functions J^nip', t) decay faster than exponentially in n: 

(1.2) \Jni^^;t)\ < Cae""" for all n. 

The need to refine this estimate will become apparent in Sect. |S1 

So far we have described asymptotic questions of a quite conventional kinship. 
The best way to introduce and motivate the new questions that we would like to 
answer, is to outline a quantum mechanical interpretation of the F-B. functions. An 
alternative physical interpretation, that considers the propagation of excitations in 
chains of classical linear oscillators, can be found in 0. 

Recall that the orthogonal polynomials {p„(/z; s)}nGN satisfy a recursion relation 
that can be written in vector form as 

(1.3) sp(Ai;s) = J^p(^;s), 

where p(/i; s) is the infinite vector of orthogonal polynomials evaluated at position 
s, and is the Jacobi matrix uniquely associated with /x (in the case when the 
moment problem is determined, of course). We can formally think of J^t as a self- 
adjoint operator acting in the space of square summable sequences, Z2(Z+) (for the 
precise treatment of this part see ^01): and consider the evolution that it generates 
via Schrodinger equation: 

(1.4) ij^m^^Mt)- 

In this equation, ip{t) is the wave-function, a vector that evolves in the space ^2(2+) 
and defines the state of the quantum system. At any time i, we can compute the 
projection of ip{t) on e„, the n-th vector of the canonical basis of /2(Z+): 



(1.5) 
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where (•, •) denotes the scalar product in Z2(Z+). 

The initial state of the evolution, ip{0), can be chosen freely. Letting it coincide 
with the first basis vector, eg, leads to the conclusion FI^ that ipn{t), the projection 
of the time evolution on the n-th basis state, can be precisely identified with J'n{p',t), 
the n-th F-B. function: 



The physical amplitudes of the quantum motion are the square moduli of the 
projections of the wave- function on the basis states of Hilbert space, |V'n(i)P- They 
are interpreted as the quantum probability to find the system in the state e„ at the 
time t. As such, they can be used to define the expected values of dynamical oper- 
ators. Unitarity of the quantum evolution operator, e~'*"'f , implies the probability 
conservation formula 



valid for all times t. This formula gives a new meaning to the analogous one already 
known for integer order Bessel functions. 

Think now of n as labelling the position in a regular one-dimensional lattice. 
Then, ■(/'(O) = eo describes a quantum system initially localized in the origin of this 
lattice, and consequently |^/'n(^)P describes the spreading of the quantum wave over 
this space. Figure 11.31 shows the initial part of the evolution in the case of the usual 
Bessel functions (for which the measure /i is absolutely continuous), and Figure ^31 
displays the same information in the case of a singular continuous Julia set measure. 
Differences between the two are apparent. 



(1.6) 



Mt) = Jn{ti;t). 



oo 



(1.7) 
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Fig. 1.3. F-B. functions \Jn{^J■',t)\^ versus time t and space n + 1, for dfi{s) 
and all three-dimensional graphs are zoomable for better viewing. 
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Fig. 1.4. F-B. functions \ J^ni^J^\ t)P versus time t and space n + 1, for a Julia set measure with 
A = 2.9. 



To gauge this spreading we utilize the moments of the position n, 

oo oo 

(1.8) Mt) := 5^n"|^„(i)l' = 



n=0 



Here, the index a takes all positive real values |12j. 

As it happens, for the singular measures that we are interested in, the asymptotic 
behavior of the position moments Va (t) is power-law, with non-trivial exponents: we 
therefore define the growth exponents (3^ (a) via the upper and lower limits 



(1.9) 



P (a) = — lim (inf) 



sup log ly^^t) 



a t^oo 



logt 



The functions (a) are also called the quantum intermittency functions. 

In the setting so defined, trivially (3^ (a) < 1, and /3*(a) = 1 in the Bessel case. 
For singular measures, bounds related to dimensional characteristics become of im- 
portance • under the sole request of existence of the orthogonal polynomials of /i, 
it is proven that P~{a) > dim^f (/i), where dim^f (/i) is the Hausdorff dimension of /x, 
and P^{a) > dimp(/i), the last quantity being the packing (or Tricot) fractal dimen- 
sion. Indeed, these theorems are even more general than required for our purpose: 
they apply to any quantum evolution in a separable Hilbert space, see the original 
references for details. 

Notice that the above bounds do not depend upon the index a. According to my 
definition, quantum intermittency is present when /3^{a) are not constant functions 
of the argument a. However strange it might seem at first, this case is typical of 
singular continuous measures supported on Cantor sets. The name of the game of 
much recent theoretical research has therefore been to study these functions, and 
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to track the origin of their behavior in the properties of the measure /z, and of its 
orthogonal polynomials. This is the problem that will be discussed in this paper. 

2. Kinematics, and expansion in orthogonal polynomials. The quanti- 
ties described in the Introduction can be obviously expressed in terms of orthogonal 
polynomials. In fact, the position moments Vait) can be written as 

oo 



(2.1) z/„(i) = ^ r^" / / d^i{s)d^i{r)e'^''-'')'pr,{^l; s)pn{fi; r) 
We are therefore confronted with the highly singular kernel 

oo 

(2.2) K^{r, s) := ^ n>„(Ai; s)p„(/i; r). 



Tl = 



When a = 0, we obtain the reproducing kernel of the orthogonal polynomials of fi: 
ifO(r,s)=<5^(r-s). ^ ^ ^ 

The behavior of individual F-B. functions can be rather erratic. The common 
procedure is then to perform a time average. Cesaro averaging is a common choice, 
but other forms of averaging work as well. For instance, Gaussian averaging, 

(2.3) Aoifm / e-l^f{s)ds, 

2iV7r J-oo 

where / is either Ua{t), or J'n{fj,;t), has the advantage of a better regularity in the 
windowing function: we have in fact 

(2.4) AGi'^a)it) ^^n" / / dfi{s)dn{r)xi/tir - s)pn{)i;s)pnilJ-;r), 

n=0 

where Xt^ i''^) — is a smooth analogue of the characteristic function of the interval 
[— For ease of notation, we use the convention u := t^^ throughout this paper. 

3. Distribution functions and lower bounds to the grov^rth exponents. 

In the study of the general problem p.9|l the consideration of a finite truncation of 
the a = moment, turns out to be useful. Define 

JV 

(3.1) VQiN.uj) ■= dfi{s)dfi{r)xu,{r - s)^ pjp;s)pn{fj,;r) . 

This is the Gaussian time average, up to time t = uj~^, of the sum of the squares 
of the first + 1 F-B. functions. Gaussian averaging is not as mandatory here 
as it is in the study of individual F-B. functions, since its regulating role can be 
also supplied by the summation over n, and yet I am not aware of any rigorous 
treatment involving only the Fourier kernel Xcj(r — s) — e~^^'^~^^/'^ . In any case, 
we shall maintain this ambiguity offering theoretical results that require averaging, 
and — at times — experimental results showing that averaging can be disposed of. 

In physical language, the discrete probability distribution |'(/;„(t)p — |i7n(/i;<)P 
(recall the normalization condition Eq. HI.?!) ) is called the wave-packet, and there- 
fore vq{N,uj) is the distribution function of the Gaussian averaged wave-packet. It 
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therefore contains all the information on this probability distribution, and a detailed 
control of this quantity, in N and lo, extends to the growth exponents. 

Typically, upper bounds on i'o(iV, have been found, yielding lower bounds on 
growth exponents for positive a. This can be easily seen by remembering that the 
quantum probability distribution |'0n(i)P is normalized by eq. (|1.7|) : limjv— >oo i^aiN, uj) — 
1 for all Lu; therefore, squeezing the head of the distribution fatten its tail. The original 
result is Guarneri's inequality f3~{a) > extended by Combes |14| to many- 

dimensional Schrodinger operators, further refined by Guarneri and Schulz-Baldes 
| 15) . and by Tcheremchantsev et al. I17j to a moment-dependent bound, in the 
form 

(3.2) p-{a)>D^,+^j-i{p). 

In the above, Dq{fi) are the generalized dimensions of the measure /i, of index q, that 
we shall define in Sect. \E\ The original hypothesis JS] that these dimensions exist 
for all g G R, and are finite for some q < 1 has been weakened to cover the 
case of the most general positive Borel measure /i. Notice finally that these bounds 
involve generalized dimensions of positive index, between zero and one. Inspection of 
the proofs reveals that this is a limitation of the technique, that deals rather crudely 
with the role of the orthogonal polynomials Pnip-] x). 




Fig. 4.1. Truncated, averaged moment vq{N,ll)) for a Julia set measure with A = 2.9. To the 
left of the figure, the region where the ansatz 14- If is well verified. The flat plateau to the right, at 
vq(^N,u)) = 1, stretches over all values of N that at time t = uj~^ have not yet been reached by the 
propagating wave. 
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4. Further lower bounds to the growth exponents. An improvement of 
these estimates is obtained if one controls the growth rate of orthogonal polynomials. 
The first attempt in this direction has been the renormalization theory of orthogonal 
polynomials of IFS measures [3 [S] that we shall meet in the following. Successively, 
the imaginative formula for the function vq^N^lo) proposed by Ketzmerick et al. |18| 
opened a different perspective: 

(4.1) vo{N,uj) N'^uj'^^'^^'\ 

where (m) is the correlation dimension of the measure /i (see Sect. El and d is 
a suitable constant that depends on properties of the orthogonal polynomials of /x, 
as we shall discuss momentarily (19) . Formula (|4.1I) can obviously be valid only for 
1 <C <C uj^^^^t^y'^. It predicts a scaling form, both in time and space, of the 
wave-packet. Pictorially, the authors of J18' say that, at fixed time, the initial part of 
the wave-packet decays with n as n'^~^. In Figure 1^1 the function i'o{N,uj) is plotted 
versus u and N, in the case of Figs. Il.ll and ll.4l We observe that the scaling (|4.1|) is 
well verified, in the region in space-time that corresponds to the decay of excitations, 
behind the wave-front. 

Starting from formula (|4.1|l Ketzmerick et al. have derived a lower bound to the 
growth exponents in the form /3(a) > D2{^j)/d. This result can be put on rigorous 
footing recalling the observation that the square moduli of all F-B. functions decay 
as p|. Therefore vq{N uj^^'^^i-''^ is a bounded function of uj. If in addition 

there exist 7 and c larger than zero such that 

(4.2) vo{N,Lo) <cN'<Lo'^^'^^'\ 



for all N and w, we can conclude that 

^2(/i) 



(4.3) (3^{a) > 



7 



for all positive values of a. In Sect. |51 we shall comment on the effectiveness of 
this bound in an exactly computable situation. Notice that our definition of 7 over- 
estimates the parameter d in Ketzmerick et al. surmise H4.1f) . and eq. (|4.2|l may be too 
crude of an estimate. Yet, in certain cases, numerical experiments as that of Fig. 14.11 
show that in a region of space-time the surmise is a good description of the function 
i'q{N,u!), and d is a good approximation of 7. As a matter of facts, the exponent 
d has a dimensional flavor, which mixes the asymptotic properties of the orthogonal 
polynomials Pnin', s) and the local properties of the measure /i. To see this, it is now 
time to briefly introduce the generalized dimensions of a measure ^. 

5. Generalized dimensions of the orthogonality measure. The spectrum 
of generalized dimensions Dq{fj,) of a positive measure /i is given, for real g ^ 1, by 
the law 

(5.1) / dfiir){fi{B^{r))y-' ~ ^ii-i)DM , 



The scaling law is made precise by taking superior and inferior limits, when w tends 
to zero, of the logarithm of the l.h.s. integral over the logarithm of lu. Of course, an 
appropriate formula exists also for q = 1. A thorough study of generalized dimensions 
is to be found in [201, EH- 
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We mention now for future reference an alternative approach to the evaluation 
of the scaling law (|5.1|l . Think of covering the support of /i by a family S of disjoint 
intervals /o-, of length l^, and measure tTo- fJ-ila)- Then, Dq{fi) is defined as the 
divergence abscissa of H{x,T,), 

(5.2) iJ(x,I]) :=^<e-«)^ 

<tGS 

when the generalized limit of finer and finer coverings is taken. 

We can now understand why the correlation dimension have a role in our problem 
|22|. Let us start from the expansion 

(5.3) AGi\JnitJ-;t)f) = J j dfj,{s)dfj,{r)xujir ~ s)pnin;s)pn{fj,;r), 

and observe that, when uj <^ ^/n, the variation of Pn{fJ-]s) over Bi^{r), the ball of 
radius uj centered at r, is negligible, so that Pn(A*; s) ~ Pni^J'', r) and 

AG{\Jni^^■,t)f) ^ jd^l{r)pl{^i■r)jd^i{s)x^{r-s) ~ j d^i{r)pl{^i-r)^i{B^{r)). 
(5.4) 

Now, the correlation dimension (m) is obtained setting g = 2 in eq. (|5.1|) . It 
so happens that the function p'^{^;r) does not alter the asymptotic behavior of the 
last integral in eq. H5.4II . and therefore -D2(a*) governs the asymptotic decay of the 
averaged square moduli of F-B. functions. Of course, this is not a substitute for 
a rigorous proof, that has been obtained in a variety of ways in the literature, as 
explained in detail in P|. 

6. Asymptotics of the orthogonal polynomials and growth exponents. 

We can now return to the wave-propagation problem, and apply the same approxi- 
mation as in eq. (|5.4|) to i'o{N,uj), to get 

N 

(6.1) MN,u;)^ d^l{r)^l{B^{r))Y,pl{^,■r). 

•' Tl = 

Suppose now that the orthogonal polynomials verify a scaling relation of the kind 

N 

(6.2) Y.Pni^^■^'■) - 9{r)N''^'-\ 

n=0 

for large N, with local dimension d{r), and a smooth function g{r) (where smooth 
is intended as a subleading behavior). Then one meets the problem, familiar in 
dimension theory, of determining the exponent d{uj), defined by 

(6.3) J dg^{r)g{r)N'^^'-^ - N'^^'^l 

in terms of the measures dg^{r) ^{B^{r))dfi{r). Suppose now that there exists 
constants C and 7 such that 

N 

(6.4) ^p2(Ai;x)<CiV^ 
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for any x in the support of then Ketzmerick et al. surmise holds. A more refined 
analysis ^1 E| can be carried on restricting the integral with respect to fi to appro- 
priate subsets fi of the support of /i, so to obtain lower boimds to 1 — vq{N,lo). This 
analysis shows that indeed under the above hypothesis 16 .41 the following lower bound 
holds: 

(6.5) /3-(a) > 

Two comments are in order: the first, is that eq. H6.5|) is a better bound than D2{fi)/'y. 
The second, that generalized dimensions of argument less than one appear. We shall 
soon return to this fact. 

In the same line is the local result [23 : suppose that there exists a Borel set S of 
positive measure, so that the restriction of to this set is a-continuous (it gives zero 
weight to any set of null a-dimensional Hausdorff measure) and so that there exists 7 
such that for any x € S H6.4|) is verified, then 

(6.6) P'{a)>-. 

1 

Further lower bounds are described in [23, [T7|. 

7. Upper bounds to the quantum intermittency function. Lower bounds 
on V[){N,Lo) do not lead to upper bounds on /3(a) [231, that arc therefore much harder 
to find [H I25L 1261 17f| , also because the strategy of restricting the consideration to a 
subset J7 of the support of /i is not sufficient here. 

The last quoted reference describes a rather interesting situation that is worth 
presenting in some detail, also because the techniques on which it is based might find 
wider applicability. One starts from the Jacobi matrices Jg.,, introduced in 28 and 
defined by 

(7.1) xpk{x;y) = {Vj^{k) + 9SQ^k)Pk{x; fi) + pk-i{x; fi) + pk+iix; ^) 

labelled by the real parameters 9 £ [6q,6i], < 60 < 9i < 00 and 77 S (0,oo). 
The structure of the recurrence relation renders Je.,, a discrete Schrodinger operator, 
with potential K,(fc). This is chosen to be null, except on a set B of selected barrier 
locations, B — {Ln,n € Z+}: V,j{k) = XB{k) k"^' ■ The exponent r] links location 
and height of the barrier. The limit rj = 00 corresponds to a Dirichlet condition at 
each L„, that clearly means no propagation and pure point spectrum. On the other 
hand, 77 = gives barriers of constant height, and generically absolutely continuous 
spectrum if these are sparse enough. Sparscness is a convenient request for analysis: 
assume that for some a > 1 and all n e Z L„+i > aL„. Under these conditions one 
can prove |27i that for all a, < a < 2, and almost all 6 one has that 

(7.2) r(«)< 



2ry + a + 1 ' 

while (a) = 1 : there is a part of the wave-packet that moves linearly in time (one 
says ballistically, in the usual jargon) while the main body follows at a slower pace. 
We refer to [2I| for the detailed analysis and illustrative pictures. 

8. Wave front propagation: unsolved asymptotics. The previous sections 
have dealt with the shape of the wave-packet behind the wave-front. This implies 
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that time, the argument of the F-B. functions, is much larger than space, the order. 
To complete the picture we must take into account what happens in the opposite 
limit, and, more importantly for our goals, in the region of the wave-front. This will 
explain our remark of Sect. O on the fact that lower bounds on vq[N,u]) in the first 
region do not yield control of /3(a). 




Fig. 8.1. Snapshots of |i/'n(tfc)P a,t exponentially spaced times t^, k = 1, . . . ,6, versus N and 
k for an I.F.S. measure, described in Sect. IISI At the bottom of the graph, the fitting line for the 
wave-front indicates the law n ~ f, with rj = 0.84165. 

A sequence of snapshots illustrating the wave-packet at exponentially spaced times 
is shown in Fig. 18.11 Over the time-span of the figure, the wave enlarges its size by 
more than two orders of magnitude. Two characteristics are to be remarked: the 
decay of the wave-packet is clearly consistent with eq. (|1.2|l . but in addition it takes 
place rather abruptly past a wave-front position. The motion of this point, on the 
other hand, appears in Fig. IS.ll to follow an algebraic law, with exponent -q. These two 
characteristics combined imply an upper bound to the growth exponents, /3^(a) < ?7, 
for all positive values of a. 

In the usual Bessel case, it is well known that rj = 1 (Bessel functions decay 
abruptly when the order exceeds the exponent), and that the wave is propagating 
linearly, so that a speed of propagation can be defined. It is interesting to remark 
that Newton's determination of the speed of sound ultimately relies on these properties 
In the singular measure case, a speed cannot be defined, and we are forced 
to introduce the intermittency function /? of the moment order a. In view of this, 
and of our previous remark on the difficulty of obtaining upper bounds to the growth 
exponents, it is of crucial relevance to develop techniques to control this particular 
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asymptotics of the F-B. functions, so to enable us, for instance, to characterize the 
exponent rj observed in Fig. 18.11 

9. Julia Set Measures: renormalization equations. We now discuss an 
example that can be worked out exactly to a large extent. We choose for /i the 
balanced measure supported on a real Julia set, generated by the quadratic map 
z ^ — A, for A > 2 The inverses of this map, 

(9.1) 0,(s)=jx/7TA 

with j = ±1, can be seen as the non-linear maps of an Iterated Function System. The 
invariant measure of this I.F.S. is defined via the equation 

(9.2) / fix)dKx) /(/ ° <l^j)i^)Mx), 

valid for any continuous function /. When A = 2, we obtain the orthogonality measure 
of the Chebyshev polynomials, suitably rescaled. When A > 2, the support of /i is a 
real Cantor set. In the graphs displayed in this paper, we have chosen A = 2.9 for no 
particular reason. 

The hierarchical structure of the support of /i is brought to evidence by iterating 
the I.F.S. maps k times: to keep the notation compact it is useful to define the index 
vector a = (cti, . . .,ak), with cr^ G { + 1,-1}, of length |(t| = fc, and the associated 
composition maps (j)cr = (t>ai ° <Pa-2 o • • • o (jj^^ . Let now /0 be the convex hull of 
the support of the measure fi, Ii^ — [—A, A], where A is the fixed point of At 
hierarchical order k = \a\, the support of fi is covered by the intervals I„, 

(9.3) h:=Mh)- 

The following remarkable property holds for the orthogonal polynomials of this 
measure [551 OT] : 

for j — ±1. Applying this property, and the balance equation H9.2|) to the Gaussian 
time averaged wave-function projections, ^G(|j7„(/i; <)p), that we denote for short 
V'^(i), we get: 



1 



(7. a' 

Iterating the renormalization procedure k times, we obtain the wave-function average 
projection at site N = n2'^', with n and k integers, in the form: 

(9.6) ?A^(t) = ^ ^ J J dfi{s)dfi{r) XujicjJair) - (j}a'{s)) PnifJ';r)pn{fJ,;s). 

(J, a' 
|<t| = \cr'\ = k 

The non-diagonal contributions, a ^ a' , in the balance equation H9.6|l have a fast 
time decay and can be neglected. In addition, in the diagonal terms, the non-linear 
maps (pa can be replaced by a linearized version: for any cr, let now 



(9.7) 



/cr(s) := 5aS + 9a 
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be the linear map, with coefficients So- and 9a-, that takes /g, the convex huU of the 
support of the measure exactly onto I^-. In other words, = 4>a{I(t>) = ^(t(-^0)j and 
the length of this cylinder is consequently proportional to 5a'- \Ia\ = "^SaA.. Usage of 
linear maps in the argument of Xuj has the effect of dividing uj hy Sa, so that 

(9.8) ^ ^'^{t6a)+£{k,n,t), 

as.t. \a\—k 

where £{k, n, t) is the error involved in the approximations we have made. The related 
error estimates are rather involved, and aim to show that £{k,n,t) is negligible, in 
appropriate asymptotic expansions. We shall boldly do this in the following. 

Eq. 1)9.8(1 is a renormalization equation, that links the wave-function projection at 
site N and time t to those at site n and earlier times. As opposed to simple estimates 
of Uo{uj,N), this equation offers us a means of controlhng the growth exponents. We 
have developed this idea in jTlHj and in the more rigorous, yet less noticed ref. [51] . 

10. Julia set measures: analysis behind the wave front. We now employ 
the renormalization analysis, eq. (|9.8|l . to compute exactly the behavior of i>q{N,uj) 
in the region behind the wave front: 

n— 1 n— 1 ^ 

z.o(A^,c^)~^2Vf2'^(a;-i)~ ^ ^ (^^c^i) = ^ -Mn,uj/Sa). 

j=0 as.t. |cr|=fc j=0 as.t. \a\=k 

(10.1) 

The form V{){N,ijj) ^ ijD2{p) gglves eq. ((10. It : upon setting iTa := 2 ^, we get 

(10.2) '^i^^ ^^'"'^''^ - 1' 

as.t. \a\—k 

that implicitly determines 7. This determination is indeed transparent: comparing 
eq. ((10.2(1 with the discrete evaluation of the generalized dimensions, eq. ((5.2(1 . we 
immediately obtain 

(10.3) D2it^) = 7 A+7(m) . 
and therefore 

(10.4) 7 = 1. 

Figure inni shows the function Bq{uj,N) := A^"iw"-°2(^)i^o(^, w), whose flat left 
piece confirms the validity of the scaling ((4.1(1 . and of the value 7=1. 

According to Sect. 01 a consequence of this calculation is a lower bound on the 
positive exponents: (3{a) > D2{p). Notice that since (m) < this bound is 

weaker that the original Guarneri's inequality f3{a) > -Di(/i) [211 ■ This fact is by no 
means accidental: information of the kind 1(4.1(1 . and more general, on the i'q(N,u;) 
for small A'' (with respect to an appropriate power of w) is not sufficient to control 
the growth exponents. 

11. Surfing the Intermittent Quantum Wave. A treatment quite analogous 
to that of the previous section can be carried out for all truncated moments of order 
a > 0: 

N „ „ AT 

Va{N,uj) := ^n"ylG(|Jn(M;On = / / dn{s) dfi{r)x:.u{r~s)Y^ n>„(^; s) p„(/^; r) . 

n=0 n=0 

(11.1) 
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Fig. 10.1. Bo(w, N) ■= N~'^uj~°2(f^)uo{N,ui) for the Julia set measure of Fig. 

Using the renormalization eq. (|9.8|l in the new situation leads to the result 

(11.2) i/„(iV,tj) - iV'^+" 

valid in the regime of decaying F-B. functions, in the leftmost part of Fig. 111.11 where 
this behavior is clearly observed. 

Notice that a new scaling region appears now to the right of the figure, ahead of 
the wave front, replacing the plateau that was obtained for a = 0: in fact, when the 
lattice site N has not yet been reached by the wave, Va{N,Lo) is independent of 
and is equal to the (Gaussian averaged) position moment of order a. Therefore, in 
such region, 

(11.3) VciN.uj) - N° uj-°'^^°'\ 

in which the intermittency function I3{a) explicitly appears. Of course, this is a trivial 
observation. It can be turned into a constructive theory only if we can stretch our 
approximations to reach this region. 

Pictorially, but appropriately, we can say that the lower bounds mentioned in 
the previous sections have been obtained by floating safely in the calm waters behind 
the wave-front. To the opposite, a complete theory of quantum intermittency can 
be obtained only if we are brave enough to catch the wave, boldly surfing on our 
approximation board the roaming waters of the F-B. wave-front, vividly depicted in 
Figure 11.41 
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Fig. 11.1. Gaussian averaged, truncated first moment ui{uj,N) for the Julia set measure with 
A = 2.9. 

Achieving this goal is a rare accompUshment: the renormahzation approach is the 
board that has enabled us to do this for Julia set measures 0i8,,34.- First, 

oo 

(11.4) u^{t)^Y.^>'{j2>^ri^%4t). 

Then, employing again the renormahzation eq. 1)9. we obtain 

C30 

(11.5) ,^a{t)^ J2 2'=^""'^E^'"^f(^-*) = 2'^""'^E^"('5-^) 

(Ts.t. \(T\—k j—0 as.t. \a\—k 

This relation has the scaling solution i^ait) ^ t^f^i*^) ^ whence by consistency 

(11.6) J2 '5-''*"^ 1' 

as.t. \a\—k 

that unveils, by comparison with eq. H5.2|l . the fundamental Julia set relation 

(11.7) /3(a) = I?i-a(M) , 

that links growth exponents and generalized dimensions. 

We have verified numerically that this relation holds exactly even without time- 
averaging [71131^. Indeed, Fig. 111.21 displays the truncated, instantaneous value of 
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Fig. 11.2. Truncated first moment vi{lo,N) for the Julia set measure with A = 2.9 without 
time-averaging. Notice the decay of fluctuations with increasing values of N . 

the first moment versus oj and N: as we have remarked previously, summation over 
n suppHes the regularizing effect. 

In a relation formally written as eq. 1)11. 7|l . but different in meaning, has 
been obtained by a renormalization procedure over Fibonacci Jacobi matrices. In 
such relation /3(a) are the growth exponents of moments averaged over initial sites 
(which mathematically amounts to averaging over different Jacobi Hamiltonians) , 
and Dq are the thermodynamical dimensions of the logarithmic potential equilibrium 
measure, that we shall also consider in the next section. A proof of this result for a 
family of Jacobi matrices has been obtained recently |87j . 

12. Linear I.F.S.: renormalization theory and a conjecture. The results 
obtained in the previous three sections are certainly neat, but by no means universal. 
They stem from the clean renormalization properties of Julia set orthogonal poly- 
nomials, eq. H9.4|l . in a situation characterized by other remarkable symmetries, the 
most notable of which is perhaps the fact that the measure of the zeros of PnifJ-', s), 
the logarithmic potential equilibrium measure, coincides with /i itself. In addition, it 
is clear that we cannot approximate an arbitrary measure with Julia set measures. 

To the contrary, linear iterated function systems "SF/SS", ID", j^J , in which we have 
at our disposal an unlimited number of maps of the kind 1)9. 7|l . li{s) = 5iS + 0i, and 
associated probabilities tt^, i = 1, . . . , Af can approximate arbitrarily well any measure 
with bounded support j42| . These I.F.S. define invariant measures ft via the obvious 
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Fig. 12.1. Distribution functions M(s) := J" d^{s') for three uniform Gibbs measures with 
the same Dq{fi) = Do = log 2/ (log 5 — log 2), as described in the text. 

generalization of eq. 19.2|l . 

r M 
(12.1) / f{s)d^l{s) = ^ TT, / (/ o k){s)d^Ji{s). 

i—1 

Moreover, linearity of the maps implies the renormalization equation 

n 

k=Q 

The coefficients F have a profound meaning, as they are the Lanczos vectors associated 
with a generahzation of the Jacobi matrix IIHI. In ref. |HI I have employed eqs. 
p2.1|l . (|12.2l) in a similar fashion than in Sect. ^] to study the intermittency function 
P{a). It has not been possible, though, to close the asymptotic relations exactly, but 
only to obtain a sequence of approximation of the intermittency function. Nonetheless, 
this approximate renormalization theory has shown that the key to the asymptotic 
behavior of the moments i^ait) lies in the properties of the coefficients F. Needless to 
say, these properties are rather elusive. 

In closing this paper I want to discuss an additional piece of evidence from 
that might give us a clue on the general problem, and is still (to my knowledge) 
unexplained. Consider the set of linear I.F.S. generated by just two maps, for which 



(12.3) 



TT, = SP>, i = 1,2 



18 



G. MANTICA 




Fig. 12.2. Intermittency function /3(a) for I.F.S. generated uniform Gibbs measures with the 
same spectrum of generalized dimensions, Dq{ii) = Do = log2/(log5 — log 2). Data plotted are for: 
two-maps I.F.S. (three lower coincident curves, crosses); a symmetrical three-maps I.F.S. (central 
curve, open squares), and an asymmetrical three-maps I.F.S. measure (top curve, filled squares). 



where Dq is a real number between zero and one, that must obviously satisfy the 
probability conservation equation 

2 2 

(12.4) l^Y.^,=Y.5f\ 

i=l 4=1 

Because of the latter equality, Dq is the box-counting dimension of the support of /i. 
Clearly, because of eqs. (|12.3|l and (|12.4|l . only one parameter among the map weights 
and contraction rates is left free, and can be put in one to one relation with Do- 

Moreover, the two affine constants Oi play no role in determining the power-law 
behavior of the moments I'ait), since we can translate and stretch linearly the support 
of fi with the only effect of multiplying the F-B. functions by a complex number of 
modulus one, and of linearly rescaling their argument |4l| . Notice finally that eq. 
()12.4|) also implies that the I.F.S. is disconnected. 

In conclusion, the family of two-maps linear IFS measures satisfying eq. \12.'J\) 
can he partitioned into equivalence classes labelled by the box-counting dimension 
D(){^). The distribution functions of three measures in the same equivalence class 
are displayed in Fig. 112.11 These measures enjoy distinctive properties. First of 
all, they are uniform Gibbs measures, according to the theory of Bowen j45j . More- 
over, since eq. ()5.2|l . with 5i and tt^ in place of la and tto-, and II(x) — 1 instead of 
the asymptotic relation, defines the generalized dimensions of linear I.F.S. measures 
exactly, one finds easily that Dq{p) = Dq{^i) for all real values of q. 



FOURIER-BESSEL FUNCTIONS OF SINGULAR CONTINUOUS MEASURES 



19 



Now, figure [TITB shows the functions (3{a) extracted numerically for three I.F.S. 
belonging to the same equivalence class Dq{p) = ^^^i^ 2 • "^^^ coincidence of the 
curves (crosses) within numerical precision — that we have also verified for other values 
of Dq — lead us to conjecture that the intermittency function /3(a) is an invariant of 
the equivalence classes defined above. 

13. Prom linear I.F.S. to potential theory: another conjecture. Even 
prior to a formal proof of the conjecture just proposed, accepting its validity leads 
to interesting speculations, and raises intriguing questions. Clearly, the conjecture 
disproves any relation of the kind (3{a) = Dq(^a)i with q a function of a, like in the 
bounds 1)3. 2|l . or in the Julia set relation 111.7|l . This is not necessarily bad news: 
it is just telling us once more that characteristics of the measure fi other than the 
generalized dimensions determine the dynamics: a few of these have been presented 
in this paper. How are the exponents d{r) and d{uj) of eqs. Hti.2|l and related, in 
and across the equivalence classes above? And the coefficients F ? Moreover, we can 
also ask how curves with different values of Dq map among themselves. But mostly, 
since eq. (|12.3|l is magnificent in its simplicity, is it there a simple argument to prove 
the conjecture? 




I I I I I I I I I I I 
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Fig. 13.1. Distribution functions M{s) := J" du(s') where v are the equilibrium measures 
associated with the three I.F.S. measures of Fig. 112.11 
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Fig. 13.2. Generalized dimensions of the three I.F.S. measures of Fig. 112. ll fsouares. horizontal 
line) and of the associated equilibrium distributions (stars) of Fio. \lS. 1\ 



En suite, notice that the extension of the conjecture to I.F.S. with three or more 
maps, without further specifications, is not vahd. In fact, an instructive counter- 
example is obtained setting M = 3, and all contraction values and weights equal 
among themselves: = 5 < i, tt^ = i, for i = 1, 2, 3. In this case, out of the 
three afhne constants 6*^, two can be set arbitrarily (for instance, so that [—1, 1] is the 
convex hull of the support of /i) , and one is left free to vary. This can be done so that 
the resulting I.F.S. is disconnected: its hierarchical structure is then composed of the 
iteration of three bands, the position of the central of which is variable. The one- 
parameter set of I.F.S. measures so obtained is composed of uniform Gibbs measures 
with Dq{fj.) = -Do(/Lt) = — log((5)/ log(3). And yet, the functions /3(a) are not invariant 
in this set: see figure ITT^ where S is chosen so to obtain the same Do as in the two- 
maps case. 

We can try an explanation of this fact. These latter three-maps I.F.S. measures 
have different quantum intcrmittcncy functions, even if they coincide "cylinderwise" , 
because their logarithmic potential equilibrium measures are different, since gaps be- 
tween covering sets la have different geometric ratios, and consequently, we expect 
different asymptotic behavior of their orthogonal polynomials. 

Having realized this, let us go back to the two-maps case. In what respect are 
then the equilibrium measures within the two-maps equivalence classes defined above 
"equivalent"? Direct inspection of their distribution functions. Fig. 112.11 provides 
no clue. It is Fig. 113.21 that contains the answer: the generalized dimensions of 
the equilibrium measures of two-maps I.F.S. in the same equivalence class are the 
same. One can therefore take the risk of putting forward a bolder conjecture: the 
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intermittency function of uniform Gibbs measures, whose equilibrium measures are 
characterized by the same spectrum of generalized dimensions, are the same |46| . 

14. Conclusions. I have discussed in this paper a number of topics that have 
originaUy been developed by mathematical physicists interested in quantum mechan- 
ics, as it appears clearly from the list of references, but that would certainly profit a 
lot from the interest of specialists in orthogonal polynomials, special functions, and 
potential theory. In fact, my formulation via Fourier-Bessel functions, the original 
idea Wl , 6 to study these problems in relation with Jacobi matrices of Iterated Func- 
tion Systems, and my introduction of the renormalization approach of orthogonal 
polynomials denote clearly how much I owe to the community that has gathered for 
this conference, and that I had the fortune to meet back in my postdoc years here in 
Atlanta. 

I have presented novel results on the asymptotic properties of F-B. functions 
for Julia set invariant measures, relating different asymptotics of the "wave-packets" 
Va{N, oj) to the properties of the invariant measure, and of its orthogonal polynomials. 
But mostly I have put forward open problems, numerical results, and conjectures that 
indicate-I hope-where to search for complete answers. How to turn this insight into 
a constructive technique for determining the intermittency function is the job that 
stays ahead. Having so arrived at the main topic of this conference, potential theory 
and its applications, I can certainly renew my best wishes to Ed, and retire in order. 

Appendix: Numerical Techniques. In addition to standard procedures, the 
research described in this paper has required novel numerical algorithms for two main 
problems: the construction of Jacobi matrices of linear I.F.S., described in ^S], |48j . 
and the computation of the F-B. functions |11) . 
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